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^2 ! Abstract 
Oh, 

O ■ Using the non-relativistic version of Borel sum rules, combined to 

. the very accurate experimental information on the Charmonium and 

Upsilon spectra, we rigorously investigate the Euclidean mass of the 
Charm and Bottom quarks. Our analysis is performed with an im- 
^ . proved expansion of the vector correlator function in the infinite heavy 

^ | quark mass limit. The optimal estimates, which take into account 

the most recent world average value of the strong coupling constant, 
as well as a conservative range of the gluon condensate values, are: 
m Eud _ 1 20 ± 0.034GeF,mf" d = 4.18 ± 0.037GeU. Their conver- 
sions to the corresponding pole mass and running mass of c and b 
quaks give respectively: M c = 1.49 ± 0.08GeV, M b = 4.65 ± OmGeV 
and m c (m c ) = 1.21 ± 0.08Gey, m b (m b) = 4.20 ± OMGeV in good 
agreement with the most recent direct estimates. 



* e- mail : ucadp @ cyb ernet.net. ma 



1 



1 Introduction: 



Heavy quark masses play a central role in the Standard Model of particle 
physics. They are often present for the phenomenology of a plethora of pro- 
cesses(weak decays, charm or bottom production,...). However, they cannot 
be estimated within the theory . In addition, because of their confinement 
at large scales, they are invisible to the experimentalists. In this context, 
different theoretical analyses have been used to the extraction of the charm 
and bottom quark masses. Besides lattice calculations, QCD sum-rules - 
a-la Shifman Vainshtein-Zakharov have proved to be an extremely useful 
framework to learn about these masses[l-5]. However as the quarks are not 
the physical states of QCD, there is no unique physical definition of quark 
mass. The most widely used definition is the pole quark mass Mq, defined 
in perturbation theory as the position of the singularity in the renormalized 
quark propagator. For the charm or bottom quark, such definition becomes 
meaningless beyond perturbation theory. Indeed it suffers from an intrinsic 
ambiguity generated by the so called infrared renormalons[17]. The latter 
induce an additional uncertainties which are presently estimated to be of the 
order of the QCD scale parameter Aqcd- 

On the other hand, the authors of ref. [1,2] observed that the radiative 
corrections are too large in the case of the heavy two point correlator associ- 
ated to the vector current Q^^Q, where Q denotes the quantum field of the 
charm or bottom quark. 

In spite of that, most QCD sum-rules determination of the c or b quark 
mass have ignored the nonperturbative renormalons effects through the use 
the pole quark mass definition. Moreover, they have performed QCD sum 
rules analysis of the heavy correlator, with the perturbative contribution 
given only up to the first order in the strong coupling constant. [] 
An alternative definition to the pole quark mass Mq is provided by the Eu- 
clidean mass rriQ Ucl , which is defined in Landau gauge at the Euclidean 
point p 2 = —Mq . Although gauge dependent, the Euclidean mass is free 
from renormalon ambiguities. It first has been introduced by SVZ [2] in 

1 Large order radiative corrections are still lacking, except a three loop contributions 
to the heavy quark correlator, calculated in part numericaly within the method of Pade 
approximants [20] . 
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order to reduce the effects of higher multi-loop corrections on heavy quarko- 
nium correlators. Using relativistic Hilbert sum-rules, they have obtained 
the following estimates: mf ud = 1.26 ± OMGeV and mf = ud 4.23 ± 0.05GeV. 
Subsequently; Reinders; Rubeinstein and Yazaki [3] confirmed the predic- 
tions of SVZ, but they have advocated the use of the Hilbert Moments at 
Q 2 = —q 2 ^ for heavy quark systems. Later on, Guberina et al. [4] and 
Reinders [5], following the same strategy, found for the b Euclidean mass: 
4.19±0.06Ge\/ and A.17±0.02GeV respectively. The discrepancies between 
these estimates come from two main sources. First, the evolution of the 
experimental information on the charmonium and the upsilon families. In- 
deed, only four resonances from six in the bottomniun were observed in the 
beginning of eighties when these papers appeared. Moreover, the measured 
hadronic parameters describing these resonances ware contaminated with 
quite large experimental errors. The second source of discrepancies is related 
to the theoretical inputs ((a s G 2 ), a Sv .) used in each work. Thus, predictions 
of mf ucl are diverse and run from 4.13 to A.28GeV, while most predictions 
of mf ucl are gathered around the value 1.26GeV. Since, on the one hand, the 
present situation of the experimental data on the J/\l/ and T systems has 
enormously improved and, on the other hand, the main input parameters of 
QCD sum-rules analysis are better under control than the past, a new and 
independent determination of the Euclidean mass of the charm and bottom 
quarks is certainly called for. 

In the following, we will present an update of the theoretical analysis of 
the vector channel of charmonium and bottomnium families based on non- 
relativistic Borel type sum-rules. We start from the expressions given in the 
original work of Bell and Bertlman [6,7], which we expand , in the infinite 
quark mass limit (rriQ Ucl — > oo), up to order. Next step will follow from 

Q 

the famous quark-hadron duality, which will permit to confront two repre- 
sentations of the ratio of Borel moments: one based on OPE calculations and 
the other on the data. An important by-product of this analysis will be the 
prediction of the Euclidean mass rriQ Ucl . 

The key input parameters most relevant to this sum-rules analysis are 
the strong coupling constant a s and the gluon condensate (a s G 2 ). For the 
former, the world average central value a s (M z ) =0.119 is by now stable and 
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the experimental error taken by all measurements ranges from 0.003 to 0.005 
[8]. In this work, our reference values will be: a s = 0.119 ± 0.005. These 
values will be combined to the formulae of the running strong coupling with 
two loop accuracy, in order to evolve a s down to the sum-rule energy scale. As 
for (a s G 2 ), first estimated by SVZ from the analysis of the charmonium sum- 
rules to be O.OSSGeV^p], there has been a lot of activity and discussion of its 
value. A number of independent estimates have appeared in the literature 
[9]. Unlike the previous quark mass determinations which use a fixed value 
of the gluon condensate, we will take the conservative range of values: 



The plan of this paper is the following. Section 2 will be devoted to review 
the most common mass definitions and their interrelations. In section 3 we 
present the quarkonium correlation function and its non-relativistic Borel 
moments. Special emphasize is put on their expansion in the infinite heavy 
quark mass limit up to Tjrorder. In section 4, rigorous analysis of sum-rules 



of the Euclidean mass of charm and bottom quarks is performed and finally, 
our numerical estimates and discussion of the uncertainties entailed by the 
sum-rules are given in section 5. 



The quark mass definition is quite confusing. Indeed, several definitions 
have been proposed and the choice of one among them is often correlated 
to the physical process of interest: each particular case requires its defini- 
tion. The most common ones rely on purely perturbative calculations. The 
pole mass Mq, the running mass in modified minimal subtraction scheme 
TnQ MS (fi) and the Euclidean mass rriQ Ucl turn out to be very popular and are 
interrelated. Indeed, it is straightforward to derive the following relations 
at- next-to-leading order in QCD coupling constant [10], 




(1) 



Q 



2 Quark mass definitions 






(3) 
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Eq. (2) relates the pole quark mass to the running quark mass defined in the 
MS scheme and normalized at the scale /x, while the connection between the 
Euclidean quark mass in Landau gauge and the pole mass is given by eq.(3). 



The running coupling constant is given by, 
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which defines the QCD scale parameter Aq C d at two loop accuracy. (3 and 
fli are the first beta functions governing the evolution of a s . They are scheme 
independent and are equal to [16], 



A) = 11 - ^n f , 



(3, = 102 - |* n/ 



(5) 



where rif is the number of flavors with mass below /j. At this stage, it 
should be understood that the number of flavors changes as \x crosses a quark 
threshold. As Aq C d depends on rif, value of Aq C d for different number of 
active flavors are defined by keeping a s continuous at the threshold scale 
/i = rriQ. Therefore, according to this prescription [11], we first extract the 

value of Aqcz) 5 '* from the present world average values of a s = 0.119 ± 0.005. 
The corresponding range is, 

Afe 5) = (235±2j)MeV (6) 



Then owing to the matching condition a s (/i = m(,) n/ _ 4 = a s (/i = rn^) n ^_ r , 

we relate Aq<£^ to Aqq^ . The corresponding range of values for rif = 4 is 
the following, 

Afe 4) = (M0±H)MeV (7) 



Equations (6) and (7) constitute important inputs in this work. 
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3 Quarkonium Sum-rules 

The basic amplitude considered in this sum-rules is the correlation function 
Ufj, v (q 2 ), induced by the electromagnetic vector current j M = Q^^Q and de- 
fined by, 

% J d 4 xexp(iqx) (SI | j M (a;)j+(0) | Q) = (q^q u - q 2 g^ Tl(q 2 ) (8) 

The structure function II (q 2 ) is related to its imaginery part through the 
dispersion relation: 

II(Q 2 ) = - f ds + subtraction, Q 2 = -q 2 (9) 

7T J (s + Q 2 ) 

where Imll(s) is proportional to the total experimental cross-section of e + e~ 
annihilition into final states with open or hidden Q flavores, 

ImIi ^ = T^ Q ^ e ^ m (10) 

cq denotes the electric charge of the heavy quark Q. In QCD sum-rules 
approach, Yi{Q 2 ) is calculated within Wilson Operator Product Expansion 
(OPE) [12], which take into account the perturbative contributions as well 
as the non perturbative effects absorbed in the vacuum matrix elements of 
operators such as G al3 G a p or qq. In the case at hand, the leading and next- 
to-leading perturbative contributions for the vector current can be found in 
Schwinger's book [15], 

ImU(s) = i-„(3 - v 2 ){l + ^[n 2 2v 2 - - 3 -)]}Q(v 2 ) (11) 



where v — (1 — 4^-)^, while the non-perturbative corrections to IT(s) will be 
represented by the leading contribution which is proportional to the gluon 
consensate (a s G 2 ) and is given by, 

n . . 1 ( 3(v 2 + l)(l-v 2 )\ .1 + v. 3v 4 -2v 2 + 3, n2 

= ^ £ ) — > < > 

(12) 
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In the sum-rules analysis, the subtracted dispersion relation in eq.(9) is im- 
proved by the Borel transform [2] leading to the exponential moments, 




ds exp(— sa)Imll(s) 



(13) 



The exponential weight in (13) has the merit to cut off large contributions and 
to take under control the effects of high dimension vacuum condensates. This 
means an improvement of the convergence of the OPE as well as a better 
enhancement of the quark-hadron duality. 

In the present work, we shall be concerned with the non-relativistic ver- 
sion of Borel sum-rules defined as, 



where we have introduced the new sum-rule parameter r = 4mQ Ucl a and the 
heavy quark energy E, through the relation s =(2mg ud + E) 2 . Moments 
present two important features: On the one hand, they can be evaluated by 
using the detailed experimental data for imIT(s)(~ cross — section). On 
the other hand, the moments in eqs. (13) and (14) are very sensitive to the 
heavy quark mass. Therefore, by equating the theoretical moments based on 
OPE calculations to the corresponding ones that use highly accurate exper- 
imental information in the heavy channel associated to the vector current , 
one can make a reliable prediction of the heavy quark mass. However, we 
will consider the ratio of the (Exponential / Borel) moments instead of the 
moments themsleves: since in contrast to the moments, the ratio has weaker 
dependance on higher order radiative corrections to the unit operator of the 
OPE . As shown in ref.[7], the non-relativistic ratio of Borel sum- rules is 
given by, 



where at r-stability region the ground state mass is, in principal provided by 



1 r 

Mr = - / 

7T JO 



dEexp(-ET)ImU(E) 



(14) 




(15) 
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the minimum of Rnr- 



The complete analytic expressions of the moments have been calculated 
by Bell and Bertlmann in terms of Whittaker functions [6]. Owing to their 
asymptotic properties, we expand them in powers of l ucl up to the next-to- 

m Q 

next-to-leading order. The non-relativistic ratio is then given by, 



/ \ c Fw/r-, 3 r , 5 10 n, Jtxuj , Nr , ,2 3 . ,107 51 . 9l 
2 In 2 , . r 5 9 , 7r 2 a s ~ 1 r 4 5 . , 

+ — a ,(.)[i + ^] + T < ^ > j^jja^P + r " j^D (16) 

where u = e uc1 . Our expression in (16) agrees with the one given by ref. [13] 

m Q t 

which take into account the mass corrections to the perturbative part of the 
ratio up to-the-first order. However, we disagree with the expansion done 
in [14]. At this stage some remarks are in order: The ratio R NR is a sign 
alternating serie, therefore, it is not always fully justified to keep only the 
next-to-leading mass corrections (NLO) and to ignore the next-to-next-to- 
leading ones (NNLO). Indeed, for comparison, we show in Table3 the results 
of the sum rules analysis for the c and b-Euclidean mass, with the ratio (16) 
expanded up to the leading order, NLO and NNLO respectively. As we see, 
in contrast to the case of the b-quark mass, the effects of the higher order 
mass corrections are important for the c-quark mass determination which is 
sensitive to the order of the l ucl expansion. 

m Q 

On the other side, the experimental ratio can be thought of as a sum 
of zero width resonances, plus a continuum contribution approximated by 
perturbation theory above a threshold St- In this approximation M(cr) is 
evaluated through, 



3 1 f 

M exp (a) = M R T R ( e+e ~ ~* R ) ex P (- aM Rl + - / ds e~ as lmllp ert (s) 

4d e 2 Q R v ' v 7r Jsy 

(17) 

where a is the QED effective coupling constant. Experimental data in J/ty 
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and T regions is very accurate. Six resonances have been observed in the 
vector channel of charmonium and bottomnium systems. Their masses and 
electronic widths are detailed in the Review of Particle Properies[18]. They 
come endowed with experimental errors (AMr, AI^) that we can safely ignor 
in our analysis. Indeed we show, for instance, that the influence due to the 
variations of Tr on the Euclidean mass estimates is very small: it amounts 
to 0.1% for the c-quark and to 0.2% for the b-quark as is indicated in Table 2. 

As for the continuum threshold, we take it somewhere above the last 
resonance of the family of interest. Then, we sensibly increase it in order 
to test the stability of our results with respect to st- For our analysis, 
the incertainties induced by the variations of st from 4.5 to 8 GeV 2 for 
charmonium and from 11.2 to 15 GeV 2 for bottomnium, are insignificant. 

4 Sum-rules analysis 

On the basis of the above discussion, we are now able to extract the Euclidean 
mass of the charm and bottom quarks by confronting the two representations 
shown in (16) and (17). Such confrontation should be realized within the con- 
text of the duality between the hadron and quark-gluon descriptions, invoked 
by SVZ. This means that the agreement between the experiment and the the- 
ory should be ensured in a compromising region of the Borel parameter r 
0, over which both contributions, perturbative and non-perturbative, to the 
theoretical ratio of moments are on the same footing. In the other hand, 
the sum-rules defined through the experimental ratio (17) should be satu- 
rated by the ground state within the Borel region: The sum-rules become 
increasingly sensitive to the ground state contribution, while the effects of 
the higher mass states are less substantial and should not exceed a given 
moderate percentage of the full contribution. 

In Figure 1, we study the behavior of both the theoretical and the experi- 
mental ratio as a function of the Borel scale a for the Charmonium family.The 
input parameters are kept fixed to their central values, < a s G 2 >= 0.55Ge^ 4 
and Aq^ 4) = 340Me1A As expected, we see that above a = Q.QhGeV 

2 This is the so called fiducial window in SVZ' works 
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the sum- rule (17) is dominated by the contribution of a single resonance 
J/ty . The analysis of theoretical expression (16) for different values of a 
(or equivalently r) shows a (Borel) region of stability where a clear mini- 
mum appears at a ~ 0.7GeV~ 2 . The best fit of (16) to the experimental 
curve is realized with m^ ucl = 1.20GeV. Next step, we study the effects 
due to the variations of the inputs on the determination of the Euclidean 
mass of Charm. With < a s G 2 >= 0.6GeV 4 , we show from Figure 3 that 
the changes in Aq^^ within the range given in (7) induce changes in mf ucl 
from 1.176 to 1.219GeV. Now, fixing the QCD scale parameter to its central 
value 340MeV, the variations of the gluon condensate in the interval 0.03— 
0.08Gel /4 generate an uncertainties on the Charm quark mass of the order 

Am Eucl = +18 Me y ( Table ^ which 

is necessary to the restoration of the 
agreement between the two descriptions. Therefore, our final results is, 

m Eucl = L2Q ± Q Q MGeV ( lg ) 

Similar analysis is performed for the Upsilon family. In figure 2, we plot the 
curve of the theoretical ratio versus the experimental data. The stability 
region shows up around a ~ 0.2 — 0.24GeV~ 2 where the best agreement 
between theory and the experiment is reached for mf ud = A.18GeV, with 
Aj^= 5) = 235MeV and < a s G 2 >= 0.55GeV 4 . The sensitivity of the b-quark 
Euclidean mass with respect to the errors of the input quantities is displayed 
in figure 4 and Table 1, from which we extract our final estimate: 

mf ucl = 4.18±0.037Ge\/ (19) 
where all the errors have been added in quadrature. 

5 Conclusions 

We have extracted the numerical values of the charm and bottom quark 
Euclidean masses from a rigorous analysis of the heavy correlator associated 
to the vector current Q'j^Q, which we combine to the precise experimental 
data on Charmonium and Upsilon systems. The analysis, performed within 
the nonrelativistic Borel sum-rules, yields to the optimal results, 
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= 1.20±0.034GeV 
= 4.18±0.037Ge\/ 



From Table 2, we show the origins of the main uncertainties which enters the 
errors on the c and b-quark masses. Besides the minor uncertainty coming 
from the experimental data on Tr, they mainly result from the variations of 
the gluon condensate in a conservative range of values and from the present 
world average QCD coupling constant at Mz scale, evolved down to the 
sum rules scale. Our predictions improve the previous determinations of 
the Euclidean mass from Hilbert Quarkonium sum-rules [2-5]. Moreover, the 
conversion of the values in (19) to the quark pole mass or MS mass of Charm 
and Beauty, through the relations (2) and (3) given to the second order in 
a s ,are in good agreement with the most recent direct estimates[19]. Indeed, 
the corresponding ranges are respectively: 



and, 



M c = 1.49 ± 0.08GeV 
M b = 4.65 ± O.OQGeV 

m c (m c ) = 1.21 ± 0.08GeV 
m b (m b) = 4.20 ± OMGeV 



(20) 
(21) 

(22) 
(23) 



where the errors on the pole masses and on MS running-masses have in- 
creased, due to the additional uncertainty related to a s (Mg) variations in 
the formulae (2,3). 
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List of Figures 



Fig. 1: Charm channel: Theoreical ratio (16) (continuous line) and the Ex- 
perimental ratio (Dashed line) as a fuction of of the Borel sum-rule parameter 
cr,with < a s G 2 >= 0.055Ge"l/ 4 , Aq^d = 235MeV and s T = A.hGeV 2 . 

Fig. 2: Similar to Fig.l in the bottom channel, with = 235MeV 

and s T = 11.2GeV 2 . 

Fig. 3: Behaviour of the Charm Euclidean mass versus the QCD scale 
parameter, wit < a s G 2 >= 0.055GeV 4 . 

Fig. 4: Similar to Fig. 3 for the bottom quark, with < a s G 2 >= 
0.055Ge\/ 4 . 
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